Frustration of tilts and A-site driven ferroelectricity in KNbOa-LiNbOa alloys 



D. 1. Bile and D. J. Singh 

Condensed Matter Sciences Division, Oak Ridge National Laboratory, 
Oak Ridge, Tennessee 37831, USA, and Department of Physics and Astronomy, 
University of Tennessee, Knoxville, Tennessee 37996, USA 

Density functional calculations for Ko.sLio.sNbOs show strong A-site driven ferroelectricity, even 
though the average tolerance factor is significantly smaller than unity and there is no stereochem- 
ically active A-site ion. This is due to the frustration of tilt instabilities by A-site disorder. There 
are very large off- centerings of the Li ions, which contribute strongly to the anisotropy between the 
tetragonal and rhombohedral ferroelectric states, yielding a tetragonal ground state even without 
strain coupling. 

PACS numbers: 77.84.Dy, 71.15.Nc 



There is current interest in discovering new Pb free 
ferroelectric and piezoelectric materials for applications. 
(K,Na)Nb03-LiTa03 [1] textured materials were re- 
cently found to have similar piezoelectric performance to 
Pb(Zr,Ti)03 (PZT)j^hich is the most used material for 
applications . [1, 0, 0] These systems are solid solutions 
with ABO2, perovskite structure similar to PZT where 
the Pb atoms at the A-site are substituted by K, Na and 
Li atoms. This raises the question if a good room temper- 
ature Pb free relaxor ferroelectric material can be found. 
Also there is theoretical interest in microscopic under- 
standing of the morphotropic phase boundary (MPB) 
between the rhombohedral and tetragonal ferroelectric 
phases, which plays a crucial role in the piezoelectric per- 
formance of materials. [HI, 0, 0] 

Lattice instabilities of perovskites are often discussed 
in terms of the tolerance factor, t = {tq + ^A)/2"'^/^(ro + 
r^), where ro, r^, and rs are the ionic radii of the O, 
A, and B ions. 0, |yj Ferroelectrics with t >1 are called 
B-site driven, because the B-site ion is too small for its 
site and can off-center. These materials are typically fer- 
roelectric with rhombohedral or similar structures and 
they do not show vertical MPB's. Examples of such ma- 
terials are BaTiOs, and KNbOs. Ferroelectrics with t <1 
are called A-site driven and they are in general not fer- 
roelectrics because of different tilts and rotations of BOq 
octahedra, which preserve the inversion symmetry. Ex- 
ceptions are materials with Pb^+ as the A-site ion, as 
in PZT, for which the ferroelectricity occurs along with 
piezoelectrically active MPB's. The lattice instabilities 
can also be discussed in terms of relative off-centering of 
A- and B-site ions in their cages. [10] This follows the 
classification in terms of t except that ferroelectrics like 
PbTiOs, with t ^1 are classified as A-site driven due to 
Pb stereochemistry. 

The piezoelectric performance depends on the ability 
of switching between rhombohedral and tetragonal phase 
under an applied external field. This ability can be con- 
trolled by substitutions of the atoms at A- and B-site 
which favors one of the phase over the other. This has 
lead to several theoretical calculations in PZT systems, 
mainly investigating the B-site alloying. |n|, 0, 0, Q 
Also recent studies of different A-site substitutions have 



showed increase in tetragonality when alloying with small 
A-site ions. [15|, IH, El, UM 

Here, we use density functional calculations to ex- 
plore the lattice instabilities of ordered supercells of 
Ko.sLio.sNbOs with perovskite structure. Since there 
are not stereochemically active ions like Pb^+ or Bi^+, 
one does not anticipate strong A-site driven ferroelec- 
tricity. However we find it due to a new mechanism- 
frustration of tilts by A-site disorder, [ill 13 There are 
not experimental results for the 50%/50% composition of 
K/Li at the A-site probably due to difficulty to synthe- 
size Ko.sLio.sNbOs in perovskite structure. (K,Li)Nb03 
around 50%/50% occupancy does not normally occur 
in perovskite structure, but rather occurs in tungsten 
bronze structure in a fairly broad composition range 
around 60% K. However, the perovskite structure is a 
higher density phase and should be accessible with appro- 
priate growth techniques, perhaps high pressure meth- 
ods. 

The calculations were performed using density func- 
tional theory within local density approximation (LDA). 
Although the energy differences discussed here are small, 
the LDA is known to correctly describe ferroelectricity 
in perovskites. Specifically, the end point compounds, 
KNbOs [19", ^0] and LiNbOg [21] are well described and 
even delicate balances like the non-ferroelectric charac- 
ter of KTaOs |23l are reproduced. We use the linearized 
augmented plane wave method (LAPW) [22] as imple- 
mented in both WIEN2K [23] and our in-house code, 
including local orbitals [24] to relax liniarization errors 
and to treat the high lying semicore states of the metal 
ions. The values of the lapw sphere radii were taken small 
for the relaxations. These values are: 2.0 a. u. for K, 1.85 
a.u. for Li and Nb and 1.5 a.u. for O. Convergence of 
the self-consistent iterations was performed for a mesh 
of (6,6,6) or equivalent grids of k points inside the Bril- 
louin zone to within 0.0001 Ry with a cutoff of -6.0 Ry 
between the valence and the core states. The number 
of plane waves used in the interstitial region is character- 
ized by a parameter RKmax=^mi^max, where R^^ is the 
O sphere radius and Kmax is the maximum plane wave 
vector. For our calculations RK^ax was taken to be 8. 
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The core states were treated fully relativistically, whereas 
the valence states were treated within a scalar relativistic 
approximation. The calculations for the Ko.sLio.sNbOa 
were done using 2x2x2 fee supercells with 10 atoms/cell 
and also using different orderings of K and Li atoms in 
2x2x2 sc supercells with 40 atoms/cell. We use the theo- 
retical lattice constant (3. 934 A) since the system consid- 
ered here has not been synthesized experimentally. 

In order to understand the lattice distortions and the 
trends in the ferroelectric instabilities we have consid- 
ered these instabilities separately in the 10 atom super- 
cell. We performed LDA calculations at the theoreti- 
cal volume for rhombohedral and tetragonal ferroelectric 
states, as well as octahedral rotations around [001] and 
[111] cubic axis. The calculations were carried out by 
imposing space group symmetries that allow the desired 
mode (centrosymmetric with 10 atoms/cell for the rota- 
tions and noncentrosymmetric with primitive cells for the 
ferroelectric states, which allow the ferroelectric displace- 
ments, but not octahedral rotations). For the primitive 
perovskite cell the octahedral rotation is a pure mode 
at the R point (^^25), while the ferroelectric modes are 
part of a manifold of transverse optic modes. We calcu- 
lated the energetics for the octahedral rotations around 
the two axes ([001] and [111]), whereas for the ferroelec- 
tric modes we have fully relaxed the atomic positions in 
the unit cells. In all the calculations, except the tetrago- 
nal ferroelectric, the lattice parameters were taken at the 
theoretical cubic values. For the tetragonal ferroelectric 
state we have varied the c/a ratio keeping the unit cell 
volume constant. 

The total energy calculations for the different in- 
stabilities give a tetragonal ferroelectric ground state 
with the two rotational distortions being strongly sup- 
pressed. Even without the strain the tetragonal fer- 
roelectric state is the ground state. This is different 
from the other known perovskite materials. The en- 
ergies relative to the undistorted cubic structure for 
the fully relaxed rhombohedral (R) and tetragonal (T) 
ferroelectric, octahedrally rotated around [111] (Rill), 
and octahedrally rotated around [001] (ROOl) structures 
are: -0.348eV/f.u., -0.408eV/f.u., -0.023eV/f.u., and - 
0.020eV/f.u. (f.u.=2Ko.5Lio.5Nb03) respectively. The 
strain calculations for the tetragonal ferroelectric state 
lower the energy even more giving the lowest value of 
-0.443eV/f.u. at 3% strain. These results are shown 
in Fig. [H The calculations show a energy difference of 
~95meV/f.u. between the strained tetragonal and rhom- 
bohedral ferroelectric states suggesting the possible pres- 
ence of an active MPB, between this composition and 
KNbOs which is rhombohedral. 

The ferroelectric distortions in the perovskite materi- 
als consist of displacements of the cations {A- and B-site 
ions) relative to the center of mass of the surrounding O 
anions. In order to see which cation displacements are 
the driving force to stabilize the tetragonal ferroelectric 
structure, we have plotted the atomic displacements of 
all the atoms relative to the O's center of mass for the 
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FIG. 1: Total energy for the strained tetragonal (T) ferro- 
electric state of Ko.sLio.sNbOs fee ordered supercell relative 
to the undistorted cubic structure as a function of the c/a 
ratio. The energy of the rhombohedral (R) ferroelectric state 
is shown in dashed line. 



lowest energy tetragonal structure (3% strained tetrago- 
nal structure) along the [001] direction. These displace- 
ments are given in Fig. [2l The Li ions have very large 
displacements. Relative to this, the Nb and K displace- 
ments are small. This is in contrast to most perovskite 
ferroelectrics, where there is significant off-centering of 
all ions. The tetragonal ferroelectric state occurs as a 
consequence of large Li displacements. This result can 
be understood in terms of frustration arising from the 
very different ionic radii of K and Li. Since the NbOe 
octahedra are relatively stiff, the coherence length of the 
octahedral rotations is much longer than that of A-site 
off- centering. [l5|, 0, [iH, [26] Because of this, the frus- 
tration of K and Li ions favors the A-site off-centering 
over the octahedral rotations. Information about the 
coherence length can be obtained from known phonon 
dispersions. In cubic KNbOs 0] sheets of instabilities 
are found whereas cubic PbTiOs [25i] has an unstable 
mode over the whole zone. In general, tilt instabilities 
exist only near the zone boundary and therefore have a 
long coherence length. Also the small values of critical 
thickness in PbTiOs thin films show the small coherence 
length of the ferroelectric state in A-site driven materi- 
als. [27] 

In order to see more clearly the individual atomic 
contributions to the tetragonal and rhomohedral ferro- 
electric states all the cations were individually displaced 
around the tetragonal [001] and rhombohedral [111] equi- 
libria while keeping all the other atoms fixed. These 
results are given in Fig. [3] and they clearly show that 
Li offcenters and favors the [001] direction over [111] 
by 73meV/(f.u.). Because of the large Li off-centering, 
Li contributes strongly to the anisotropy between the 
tetragonal and rhombohedral ground states, yielding a 
tetragonal ground state. This is in contrast to the usual 
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FIG. 2: Atom displacements relative to the center of mass of 
the O atoms for the 3% strained tetragonal (T) ferroelectric 
state of Ko.sLio.sNbOs fee ordered supercell. 




FIG. 3: Energy vs individual atomic displacements around 
T [001] and R [111] directions of Ko.sLio.sNbOs fee ordered 
supercell. The values are shown in: o for Li along T, □ for Li 
along R, for both Nb atoms along T, y for both Nb atoms 
along R, A for one Nb atom along R, + for K along T, and ★ 
for K along R directions. 



case in perovskite ferroelectrics where strain coupling 
may stabilize the tetragonal state, but otherwise the 
anisotropy comes from the B-site and favors rhombohe- 
dral states. This is significant because having a material 
with a true tetragonal state is important for obtaining 
a vertical MPB when alloying with a rhombohedral ma- 
terial, as is the case in PZT. Individual displacements 
of Nb and K atoms do not show any instability toward 
ferroelectric states. 

To see the effect of different K and Li arrange- 
ments on the ferroelectric state we performed full re- 
laxations on different K/Li ordered 2x2x2 supercells 
with 40 atoms/cell. We have considered four different 
K/Li ordered supercells : fee order, chain order along 
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FIG. 4: Unit cells for the 40 atom supercells with different 
orders: (a) fee order, (b) chain order along z-axis, (c) layer 
order perpendicular to the z-axis, and (d) the most disordered 
case for this type of supercell. The K, Li and Nb ions are 
shown in large gray, small gray, and black colors, respectively. 
For the reason of clarity O anions are not shown. 



the z-axis, layered perpendicular to the z-axis, and a 
disordered case for this type of supercell. These are 
shown in Fig. [H The total energies of the fully re- 
laxed cells relative to that of the supercell with fee or- 
der are: 0meV/(f.u.) for fee order, -107.12meV/(f.u.) 
for chain order, +88.54meV/(f.u.) for layer order, and 
-17.7meV/(f.u.) for the most disordered case. The most 
stable ordered configuration is the chain order. However, 
if this material follows the pattern of other perovskite 
oxides, samples will be disordered as the cations become 
mobile only at high temperature. The supercells give 
ferroelectric states and not octahedral rotations as the 
ground states with Li ions having the largest displace- 
ments. For the fee order and chain order supercells, Li 
displacements are along the tetragonal [001] direction. 
For the chain order the Li displacements (~1.5A) are 
larger that those of the fee order, whereas the K displace- 
ments are smaller. The chain order better minimizes the 
strain by larger Li displacements yielding the lowest en- 
ergy. 

The most disordered case (Fig. |3]i) for which the 3D 
plots of the Li displacements are shown in Fig [5] is partic- 
ularly interesting. Three of the Li ions off-center almost 
along the three tetragonal directions and one Li ion off- 
centers to a more general direction given by the displace- 
ment vector [0.801,0.380,0.467] (Fig [5]). The direction of 
movement for the three Li which displace close to the 
tetragonal directions is toward the next nearest neighbor 
Li ions. The Li displacements show frustration of the Li 
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FIG. 5: Li displacements relative to the center of mass of the O 
atoms for the disordered supercell. 

in finding the optimal direction of off-centering. Such a 
property is favorable for forming a relaxor ferroelectric. 

To summarize, lattice instabilities of supercells of 
Ko.sLio.sNbOs were investigated using electronic struc- 
ture calculations. Frustration between large K and small 
Li atoms at the A-site leads to strong suppression of the 
rotational instabilities and the presence of A-site driven 
ferroelectricity with large Li off-centering. This favors 
the tetragonal ferroelectric ground state even without 



strain. This is different than the other perovskite materi- 
als which have a rhomohedral ferroelectric state without 
strain inclusion. 

A-site driven ferroelectrics like PZT also have ferro- 
electric coherence lengths shorter than the B-site driven 
materials. ^ This and the coupling of Li off-centering 
to the local Li-K ordering suggest that materials near 
this composition may be good relaxor ferroelectrics. 

Ko.sLio.sNbOs system is an example of alloying large 
and small A-site cations in order to create strong A- 
site ferroelectricity with a tetragonal ferroelectric ground 
state. This is significant because it suggests ways to find 
Pb free perovskites with nearly vertical MPB's. 
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